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Abstract. We study the rate of convergence for (variational) eigenvalues of 
several non-linear problems involving oscillating weights and subject to differ- 
ent kinds of boundary conditions in bounded domains. 



1. Introduction 

In this work we study the asymptotic behavior as e — > of the (variational) 
eigenvalues of 

(1.1) -div(a(x,Vu e )) + V{^)\u £ \ p - 2 u e = A £ p(f)|u £ | p ~V in CI 

with different boundary conditions (Dirichlct, Neumann, etc.), where f2 C is a 
bounded domain, e is a positive real number, and \ e is the eigenvalue parameter. 
We also consider eigenvalue dependent boundary conditions, 

-div(a(x, Vu £ )) + V(f)K| p_2 « s = A e p(§ )\u E \ p ~ 2 u e in 0, 
a(x, Vu e )v = \ e \u e \ p ~ 2 u e in dCl. 

and the Steklov problem 

-div{a{x, Vu £ )) + V(f )\u £ \ p ~ 2 u e = in ft, 
a(x, Vu e )v — \ E \u e \ p ~ 2 u £ in dfl. 

The weight function p(x) is assumed to be bounded away from zero and infin- 
ity, the potential function V{x) is bounded and the operator a{x, £) has precise 
hypotheses that are stated below, but the prototypical example is 

(1.2) a(x,£) =A(x)\£\ p - 2 £, xefl, £eW N 

with 1 < p < +oo, and A(x) is a uniformly elliptic matrix. 

The problem of finding the asymptotic behavior of the eigenvalues of (|l.ip has 
relevance in different fields of applications, and it is an important part of Homoge- 
nization Theory. 

We will consider the important case of periodic homogenization, i.e., the func- 
tions p(x) and V(x) are assume to be Q— periodic functions, Q being the unit cube 
in R N . 

The natural limit problem of (jl.lj) as e — > is given by 

(1.3) - div(a(x,\7u)) + V\u\ p - 2 u = \p\u\ p - 2 u in Q 
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with the corresponding boundary condition, where p and V are the averages of p 
and V in Q, respectively. 

In this work, we focus our attention on the order of convergence of the eigenval- 
ues, that is, an estimate of |A| — in terms of k and e, where A| and A& are the 
A:-th variational eigenvalues of problems and (|1.3[) respectively. 

The homogenization problem for eigenvalues has deserved a great deal of atten- 
tion in the past, specially in the linear case, that is, problem (jl.ljl with a(x, £) given 
by (O and p = 2. 

In particular, the first result on the order of convergence for the linear problem 
complemented with homogeneous Dirichlet boundary conditions can be found in 
[17] where it is proved that 

|A% — A fc | <Ckv £ ^. 

where C is a positive constant independent of k and e 

Later on, again in the linear case and with homogeneous Dirichlet boundary 
conditions, Santonsa and Vogelius [21] proved, by using asymptotic expansions, 
that 

\K~M <Ce 

where C depends on k. 

We want to remark that in the above mentioned works, the authors allowed for 
an e dependance on the diffusion matrix in (jl.2l) . 

More recently, the linear problem in dimension N = 1, and for the operator 
a(x,t;) = £, was studied by Castro and Zuazua in [SJ [5J. In those articles the 
authors, using the so-called WKB method which relays on asymptotic expansions 
of the solutions of the problem, and the explicit knowledge of the eigenfunctions 
and eigenvalues of the constant coefficient limit problem, proved the bound 

|A|-A fc | < C7e 4 £. 

Let us mention that their method needs higher regularity on the weight p, which 
must belong at least to C 2 . 

Recently, Kenig, Lin and Shen |14| studied the linear problem in any dimension 
(allowing an s dependance in the diffusion matrix of the elliptic operator) and 
proved that for Lipschitz domains SI one has 

|A% — A*| <Cs\\oge\i+° 

for any a > 0, C depending on k and a. 

Also, the authors show that if the domain O is more regular (C 1,1 is enough) 
they can get rid of the logarithmic term in the above estimate. However, no explicit 
dependance of C on k is obtained in that work. 

In the non-linear case without dependence on e in a(x,£), and N > 1, we proved 
in [12] (by means of a precise order of convergence for oscillating integrals), that 

|A|-A fc | <Ck%e 

with C independent of k and e. 

Moreover, in |11| . for the one dimensional problem, we show that the constant 
entering in the above estimate can be found explicitly and, moreover, a dependence 
on e on the operator a(x, £) was treated. 

Let us stress the fact that all the above mentioned works deal with the homoge- 
nous Dirichlet boundary condition case, with the exception of the aforementioned 
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paper |14| where also it is considered Neumann and Steklov boundary conditions 
and similar results as in the Dirichlet boundary condition case were found. 

Some results are known in the linear case when different boundary conditions are 
considered. In the one-dimensional case with e dependence in the operator a(x, £) 
Moskov and Vogelius [16] by using the Osborn's eigenvalues estimate proved that 

\X%-X k \<Ce 

where C depends on k. 

Recently, in by similar methods to those in [T^] the rate of convergence of 
the first non-trivial curve in a weighted Fucik problem with Neumann boundary 
conditions was found. As corollary it follows that 

\\\-\i\<Ce, |A| - A 2 | < Ce 

where Af and A| are the first and second eigenvalues of the nonlinear equation 
with Neumann boundary conditions and C does not depends on e. 

There is a large class of usually studied eigenvalue problems related to problem 
with different boundary conditions. In this paper we analyze the most common 

ones. 

Let us consider p, V two Q— periodic functions, being Q the unit cube in . 
which satisfy 

(1.4) < p~ < p(x) <p+<+oo a.e. fi and V G L°°(f2). 

for certain constants p~ < p + . We denote p E (x) := p(f ) and V e (x) := V(|). 
For each e > fixed, we define the following eigenvalue problems: 

• Dirichlet problem: 

' -div(a(x,\7u E )) + V s \u E \P- 2 u e = \ e p E \u E \P- 2 u E in fl, 
^u £ = in an. 

Neumann problem 

'-div(a(x, Vu E )) + V e \u e \ p - 2 u e = A e p £ \u e \P~ 2 u e in 0, 
a(x,Vu e )v = in 9f2. 

• Robin problem 

' -div(a(x,Vu E )) + V E \u e \P- 2 u £ = X E p £ \u E \P- 2 u e in Q, 
a{x,Vu £ )v + l3\u e \P- 2 u e = in dil. 

• Non-flux problem: 

!-div(a(x,\7u 6 )) + V E \u E \P~ 2 u e = X e p £ \u e \P- 2 u E in 0, 
u e = constant in <9f2, 

J d(l a(x, Vw E )i/ dS = 0. 

• Eigenvalue dependent boundary condition 

' -div{a{x, Vu e )) + V E \u £ \P' 2 u e = X e p £ \u e \ p - 2 u E in O, 
a(x, Vu E )v = X e \u e \p- 2 u e in dfl. 



(1.5) D e (Q) 



(1.6) N £ (n) 



(i.7) R s (n) 



(i.9) s e (r>) 
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Steklov problem 

' -div(a(x,\7u £ )) + V E \u E \ p - 2 u e = in 0, 
a{x, Vu £ )f = X E \u £ \ p - 2 u E in 



(i.io) s e (n) 



We consider to be a bounded domain in M. N , the smoothness will be precised in 
each case. The operator a(x, £) satisfy properties (H0)-(H8) given in §2. Observe 
that in problems (|1.5[) - (|1.10[) the operator a(x,£) does not depends on e, the de- 
pendence on the parameter appears only in the weights and potential functions. By 
v we denote the outer unit normal vector with respect to dCl. The parameter (3 in 
the Robin problem is in [0, oo). We observe that when (3 = it corresponds to the 
Neumann problem and when (3 = oo it corresponds with the Dirichlet problem. In 
the boundary condition problem (|1.9p and in the Steklov problem (|1.10[) we require 
the potential function to be strictly positive, i.e., there exists V~ > such that 
V~ < V(x) a.e. in ft. Observe that this requirement is not necessary in (jl.5p ~ p.8p 
since the hypothesis of V being bounded below away from zero can be assumed 
without loss of generality. 

In (jl.5P (jl.lOp the natural limit problems (as e — >• 0) are the analogous ones with 
the weights p e and the potentials V e replaced by their averages in the unit cube Q, 
i.e. 



P= P(y)dy, V= V(y)dy. 

JQ JQ 

It is not difficult to see, for any of the problems (|1.5p - (|1.10p . that if A £ is a con- 
vergent sequence of eigenvalues as e — >• then A = lim e ^o A e is an eigenvalue of the 
corresponding limit problem and, up to some subsequence, the associated eigen- 
functions u 6 converge weakly to an associated eigenfunction u of the corresponding 
limit problem. 

For the Dirichlet problem (|1.5p this fact was proved in [T| (see also [12] for a 
simplified proof of this result). The proofs for the others problems (|1.6p - (|1.10p are 
analogous. 

Our aim is to study the order of convergence of the eigenvalues of problems 
(|1.5|) - (|1.10|) to those of the limit equations. 

Using results concerning to oscillating integrals, we prove our main results: 

Theorem 1.1. Let A| be the k-th variational eigenvalue associated to any of the 
problems (jl.5p - (|1.9p . respectively. Let A^ be the k-th variational eigenvalue as- 
sociated to the correspondent limit problem. Then there exists a constant C > 
independent of the parameters e and k such that 

\X%-X k \<Ck^£. 

Theorem 1.2. Let A| be the k-th variational eigenvalue associated to equation 
(jl.lOp . Let Afc be the k-th variational eigenvalue associated to the correspondent 
limit problem. Then there exists a constant C > independent of the parameters e 
and k such that 

\Xl-X k \<Ck^Te. 



Remark 1.3. Let us note that, for problem (|1.9p . the bound can be improved when 
p < N, and we get 



2(p-l) 

|A|-A fc | <Ck~^e. 
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The rest of the paper is organized as follows. In Section §2 we introduce the 
class of operators considered and the hypotheses on the functions a(x, £), p, V, and 
the notation which will be used. In Section §3 we analyze the eigenvalue problems 
and we show the relationships between them. Section §4 is devoted to oscillatory 
integrals, and in Section §5 we prove the main results. 

2. Preliminary results 

2.1. Monotone operators. We start this section by making the precise assump- 
tions on the function a(x, £) 

Let C R N , N > 1 be a bounded domain. We consider a : CI x l w -> R N that 
satisfies the following conditions: 

(HO) measur •ability: a(-, •) is a Caratheodory function, i.e. a(x, ■) is continuous 

a.e. x £ CI, and a(-,£) is measurable for every £ £ R w . 
(HI) monotonicity: < (a(x, £i) — a(cc, ^))(^i — £2)- 
(H2) coercivity: a|£| p < a(x, 
(H3) continuity: a(x,£) < 

(H4) p— homogeneity: a(x,t£) = t p ~ 1 a{x,£) for every £ > 0. 
(H5) oddness: a(x, — £) = — a(x, £). 

Let us introduce ^(x, £1, £2) = a(x,£i)£i + a(x,£.2)£,2 for all £1,^2 £ R^, and all 
x £ CI; and let 5 = min{p/2, (p — 1)}. 

(H6) equi- continuity: 

\a(x, a) - a(x, 6)| < c*(x, fc)^"^^, 6) - a(x, - 

(H7) cyclical monotonicity: ^i=i a ( x ; — £i) < 0) f° r au & > T an d 

Ci, ■ • - ,4+1) with 6 = £fc+l- 
(H8) strict monotonicity: let 7 = max(2,p), then 

<*|£i - 6P*(z, 6, 6) 1_(7/p) < Hx, 6) - a(x, - 

Hypotheses (H1)-(H3) are necessary to ensure the G— convergence of the op- 
erators associated to a(x,£). On the other hand, hypotheses (H4)-(H7) are all 
important in the context of a well-posed eigenvalue problem. We assume (H8) for 
technical reasons. 

We add that the conditions (H0)-(H8) are not completely independent of each 
other. It can be seen easily that (H8) implies (H1)-(H2) and that (H4) implies 
(H3) in addition to the continuity of the coefficient, for details see pQ. 

Remark 2.1. The prototype for such functions is a(x,£) = A(x)\£] p ~ 2 ^, where A(-) 
is a measurable function with values in the set of N x N symmetric matrices which 
satisfies 

a'lCI 2 <A(x)£-£, \A(x)(\< V£ £ R N , a.e. x £ Cl. 

for some positive constants a' and /?'. 

In particular, under these conditions, we have the following Proposition due to 
Baffico, Conca and Rajesh [I] 
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Proposition 2.2. Given a(x,£) satisfying conditions (H0)-(H8), there exists a 
unique Caratheodory function <!> which is even, p— homogeneous strictly convex and 
differ entiahle in the variable £ satisfying 

(2.1) <$(x,0<m P 
for all £ £ M. N a.e. x £ f2 such that 

V f $(z,£) =pa(x,0 
and normalized such that $(x, 0) = 0. 

Proof. See Lemma 3.3 in [T]. □ 

2.2. General hypotheses and notation. Throughout the paper the following 
hypotheses and notation will be used: 

• By Q we always mean the unit cube in R w , i.e. Q — [0, l] N . 

• The functions g,p, V will always refer to functions in L°°(R N ) that are 
Q— periodic. 

• For e > we denote 

«fe(aO = fl(f). Pc(x)=p(%), V e (x) = V(£). 

• The average on Q will be denoted by 

9=1 g(y) dy, p= p(y) dy, V= V(y) dy. 

JQ JQ JQ 

• Observe that g e — 1 g, p £ p and V £ — ^ V (as e — > 0) weakly * in L 00 . 

• On the weight function p we assume that it is bounded away from zero. 
That is, there exist constants < p~ < p + < oo such that 

P~ < P(x) < p + - 

• On the potential function V, for the problems (|1.9p and (jl.lOp . wc assume 
that the potential function is strictly positive, i.e., there exists V~ > such 
that 

V~ < V{x) a.e. in O. 

• Lipschitz regularity of the domain il is assumed in problems (|1.5[) and (|1.8D , 
and C 1 regularity in all other cases. 



3. Eigenvalues 

In this section we define the variational eigenvalues of problems (|1.5j) - (|1.10D 
and collect some properties and relations between them. We focus only on the 
properties of the eigenvalues that will be used later, a more detailed study can 
be found in the following references: the Dirichlet problem (|1.5p was studied by 
Garcia- Azorero and Peral Alonso in [5]; the Neumann problem can be found in 
the work of Huang [TUj; the Robin problem (|1.7[) can be found in [TS], together 
with the Non-flux problem (|1.8[) ; both generalize periodic and separated boundary 
conditions in classical Sturm Liouvillc problems, studied by several authors, see 
the paper of Binding and Rynne [3] among others; the problem with Eigenvalue 
dependent boundary conditions (|1.9[) was studied by Binding, Browne and Watson 
[2] in the one dimensional case, and the Steklov problem (|1.10j) was considered by 
Fernandez Bonder and Rossi in [4] . Finally, more general monotone operators than 
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the p— Laplacian, like the ones we will consider here, were studied by Kawohl, Lucia 
and Prashanth, see [13]. 

First of all, we observe that by replacing A e by A e + \\V\loo + V- in ([T3|h([L"g|) 
we can assume that the potential function verifies that 

V[x) > V- > 

In problems l|1.9|) - (|1.10|l this has to be imposed on V. 

By means of the Ljusternik-Schnirelmann theory (see (52] for instance) we know 
that the variational spectrum of these problems consists in countable sequences of 
positive eigenvalues tending to +oo. Define the followings functionals 



F(u,p)= [ P \u\p, 
Jo, 

G{u) = I $(i,V«), 



<2 



H(u) = / \u\p, 

where Q(x,t;) is the potential function given in Proposition 12.21 Then, for each 
fixed e > we can give the characterization of the fc-th variational eigenvalues of 
(Oll- lHTnil as follow: 



AjP = inf sup 



G(u)+F(u,V) 



(3.1) 



cef k uac F(u, p) 

N . G(u) + F{ Ul V) 
At, = mf sup —. r , 

cer kuec F{u,p) 

inf mu) + G(u)+F(u,V) 
cer kuec F{u,p) 
G{u) + F{u,V) 



AjT = inf sup 



ceT kueC F(u,p) 

B . G(u) + F(u, V) 
X? = mf sup — — ; fi 

k Cer kue % H(u)+F(u,p)' 

xS . . G(u) + F(u, V) 
\i = mf sup— — — — ^ — -. 

c^T ku £ H{u) 



Here, 



T fc = {C C W hp (n) : C compact, C = -C, 7 (C) > k}, 

f fc = {C C W ( ] tP (n) : C compact, C = -C, 7 (C) > k}, 

f fe = {C C Wq' p (Q) © R : C compact, C = -C, 7 (C) > k} 

and 7(C) is the Kranoselskii genus (see for instance [7] for definition and prop- 
erties). 

From the variational characterization, we immediately obtain the following in- 
equalities 

(3.2) Af < Af < min{Af , A*} < max{Af , A*} < Af 

(3.3) Af < Af 
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In the followings Lemmas we give upper bounds for the eigenvalues Xk defined in 
(|3.1|) in terms of k and f2. Here and in all the paper we will consider that fi C 
is a bounded domain. These estimates will be useful to prove the main results since 
it provide us with the growth rate of the eigenvalues. 

Lemma 3.1. Let Xf? , Ajfj, AjP, X^ and X^ be the k-th variational eigenvalues 
defined in ()3.ip . Then 

Af < A? < min{Af , A«} < max{Af , A«} < Af < Ck^ N 
where C depends only on £1 and the bounds (|1.4[) . (|2.1[) . 



Proof. From (|3.2j) . it is enough to prove the last inequality. 
Now, from (|2.1j) we have 

G{u) + F(u, V) max{/3, V+) J n |Vu|* + |tt|* 



F(u,p) p- J n \u\v ' 

from where it follows that 

D max{/3,y+} 

^fe S = Pk, 

P 

where pk is the k— th eigenvalue of 

j-A p u+ |u|p- 2 u = p\u\P- 2 u inQ 
1 u = on <9f2. 

Observe that it £ W ' p (fi) is solution of (|3.4p if and only if u is solution of 

{—ApU = p\u\ p ~ 2 u in O 
it = on dfl, 

where p = p — 1, which satisfies that (see [9]) 

(3.5) p k < Ck?/ N , 

and the result follows. □ 

Lemma 3.2. Let X^ and A^ be the k-th variational eigenvalue of B{VL) and S(il) 
respectively. Then 

Af <Af <Ck^h, 
where C is a constant depending on V + , a and fl. 



Proof. From (|2.1[) and (|1 .4|) we have 



G(u) + F(u,V) ^_ rl „ +1 /„|Vu|f + / n 



<max{i 



from where it follows that 

(3.6) Af <max{i V + }p k , 

where pk is the k— th eigenvalue of 

{-A, D u + \u\p- 2 u = inf2 



\\J U \P' 2 ^ = p\u\P~ 2 u in dfl. 
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Ii is proved in [18] the following estimate for \ik 




where c is a positive constant independent of k. 



Finally, from (pHl) . (jcil))) and (|3~7)l the result follows. 



□ 




4. Preliminaries on oscillatory integrals. 



In order to deal with the rate of convergence of the eigenvalues, the main tool 
that we use is the study of oscillating integrals. These will allow us to replace an 
integral involving a rapidly oscillating function with one that involves its average 
in the unit cube. 

Let p be a Q-pcriodic weight. It is well-known that p(f ) converges weakly* in 
L°° to its average over Q. We are interested in the rate of the convergence in terms 
of e. In [TJ] it is proved that 

Theorem 4.1 ([IS], Theorem 5.5). Let f2 C be a bounded domain with Lipschitz 
boundary and let g £ L°°(M. N ) be a Q— periodic function, Q being the unit cube in 
WL N . Then, there exists a constant C depending only on p, fl and H^Hl^r") such 



With the aid of Theorem 14.11 in [12] we were able to analyze the Dirichlet 
boundary condition case. To deal with different boundary conditions, we need a 
similar Theorem that allows us to include the function space W 1,p (f2). 

The fact of enlarge the set of test functions is reflected in the need for more 
regularity on the domain £1. In |20] the following result is proved. 

Theorem 4.2 ([20], Theorem 4.3). Let Q C W l be a bounded domain with C 1 
boundary and let g £ L°°(M. N ) be a Q— periodic function, Q being the unit cube in 
K . Then for every u £ W 1,p (£l) there exists a constant C depending only on p, f2 
and ||<7||l~(r«) such that 



Now, we need a couple of simple technical lemmas which are used in the proof 
of our main results. 



that 




for every u £ W ' p (Ci), where g is the average of g over Q. 




Lemma 4.3. Let u £ W X ' P (VL). Then 

w\un w ^ (n) <p\\u\\v 



W 1 .p(f2) 



<C(G(u)+F(u,p)) 
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where p is an arbitrary weight satisfying (|1.4[) and C is a constant depending on 
p, p and the bounds of (|2.1[) . 



Proof. By using Young's inequality 

IIH , 'lkM(n) = IIH , 'll^(n)+p|||«r 1 V«|| i i (n) 
(4.1) <p\\u\\l P(n) + \\Vu\\l Pm 

<p\W\\w^v { ny 

Moreover, by (|2TT|) 

\\u\\ p wl , pm <Up [ |u|* + £ / $(x,Vu) 
^ 4 2) P \ Jn a Jn 

< imax{|,l}(G(u) + F(u,p)). 
P 

From (|4.ip and (|4.2|) . the Lemma follows. □ 

Lemma 4.4. Let u G W 1,p (f2) and p, V € L°°(R Ar ) oe two Q— periodic functions 
satisfying (|1.4[) . TTien t/iere exists a constant c depending only on p, \\p\\l°°(vl n )> 
\l=°(Vl n ), the constants in (|1.4[) and (|2.1[) smc/i i/iai 

F(n,p) F(u,V)+G(tt) 
< 1 + ce- 



F{u lPe ) ~ F{u,p) 
and 

F{u lPe ) <1 + c£ F(u,V E )+G(u) 



F{u,p) ~ F(u,p e ) 
Proof. Applying Theorem 14.21 we obtain that 

(4.3) ^< 1 + CE «±a. 

By Lemma 14.31 and (|1.4j) we bound (|4.3j) as 

(4.4) 1 + C£± F(u,V) + G(u)^ 

p F{u,p) 
Similarly, by Lemma |4~31 and (|1.4[) we get 

<l + g£ max { l,n^- G(u) : FKK) . 

p 7" F(u,p e ) 

This completes the proof of the Lemma. □ 

5. Main results. 

The proofs of Theorems 11.11 and 11.21 follow the same general lines of the ones of 
Theorem 5.6 in [T^]. The fundamental tool to estimate the rates of convergence of 
the eigenvalues is the error bound of oscillating integrals given in Theorems 14.11 and 
14.21 Observe that the regularity of the domain 51 considered in equations (|1.5|) 
(|1.10|) are the necessary to apply Theorems 14.21 and 14. 11 i.e., Lipschitz regularity in 
equation (|1.5|) and C 1 regularity in all other cases. 
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5.1. Proof ofTheorem ll.il The Dirichlct boundary condition case, i.e. problem 
D £ (fl), was treated in [12l Theorem 5.6]. We prove the result in detail for problem 
N e (Q). For R £ (il), P e (f2)and B £ (£l) the proofs are very similar and we will make 
a sketch highlighting only the differences. 

For simplicity wc will denote A| and X k (without the superindex N) to the /c-th 
variational eigenvalue of N e (Sl) and its limit problem obtained as e — > 0. 

Let 5 > and let Gg C W 1,p (fi) be a compact, symmetric set of genus k such 
that 

(*u ^ G(u)+F(u,V) 

5.1) X k = sup — — =r h 0(5). 

«gg* F (u, P) 

We use now the set Gg, which is admissible in the variational characterization 
of the fc-th eigenvalue A|, in order to find a bound for it as follows, 

(5.2) A| < sup ^)+ nu V,)F(ur P ) 

y ' k ~ueGl F(u,p) F(u, Pe ) 

Now, we look for bounds of the two quotients in (|5.2p . 

For every function u £ Gg C W 1,p (Vl) we can apply Theorem 14.21 and we obtain 
that 

G(u) + F(u,V £ ) G(u) + F(u,V) \\\u\n w ^ m 
{ ' F(u,p) ~ F(u,p) +L£ F(u,p) ■ 

By using Lemma 1431 we have for each u £ Gg there exists some constant c > 
such that 

\\\u\ p \\w^(n) < G{u) + F(u,V) 
F(u,p) ~ F(u,p) 

(5.4) G(«) + F(t,, 

5= o sup _ 

= C(A fc + 0(<S)). 

Since ueGjc M /1,p (fi), by applying Lemma [4T4l and (|5.4p we obtain that 

F(u,p) <1 + c JM P Ww^u) 



(5.5) F(u,p e ) ~ F(u,p) 

<l + Ce(\ k + 0(5)). 

Then, combining (|5.3|) . (|5.4[) and (|5.5[) we find that 

A| < (A* + 0(5) + Ce(\ k + 0(5))) (1 + Ce(X k + 0(5))) 
Letting 5 — > we get 

(5.6) A| - A fc < Ce(\ 2 k + X k ). 

In a similar way, interchanging the roles of X k and A|, we obtain 

(5.7) A fc -A| <Ce((Xl) 2 + Xl). 
So, from (|5.6p and (|5.7[) . we arrive at 

|A| - A fe | < Cemax^ + A fe , (A|,) 2 + A|}. 
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In order to complete the proof of the Theorem, we need an estimate on Afc and 
A|. By Lemma 13. II we can compare them with the fc-th variational eigenvalue of 
the p— Laplacian obtaining 

\\%-\ k \<Cek 2 ^ N , 

and the proof is complete. □ 

In the remaining of this subsection, we highlight the difference between the 
Neumann case and the rest of the boundary conditions with the exception of the 
Steklov problem that has a separate treatment. 

Dirichlet: As we mentioned before, this problem was addressed in (T2J Theorem 
5.6]. Here, functions are taken in W ' p (f2) instead iy 1,p (f2) in the variational char- 
acterization of the eigenvalues. This leads to use Theorem 14.11 instead of Theorem 
14.21 to estimate the oscillating integrals. Now, for each function u G W^' p (fi) we 
can apply Theorem 14.11 and obtain an analogous equation to (|5 .3[) 



G(u)+F(u,p e ) G(u) + F(u,V) H v "IIlp ( o) 
F{u,p) ~ F(u,p) £ F(u,p) • 

Now, the difference with the Neumann case is the way we bound the quotients 
\\Vu\\ p LP{n) /F(u,p) and || Vu\\ p LP[n) / F{ Ul p s ). 
By {ft) and J2T) we get 

l|V-l| P LP( o) K P HVu||£ p(n) 



F(u,p e ) p F{u,p) 

p 1 G(u)+F(u,V) 



< 



and 



p a F(u,p) 



|V^||£ P(n) p + \\Vu\\l P(n) 



F(u,p) p F(u,p e ) 
Taking into account this changes, the proof is analogous to the Neumann one. 

Non-Flux: Let u G W 1,p (fi) © R, the n u = v + c where v G W ,p (fl) and c 
is a constant depending on u. It follows that u — c G Wo' p (n). Observe that if 
u G W 1,p (fi) © R then \u\ p G Wq'^O) © R and then |u| p - c G Wq'^O). By using 
Theorem 14.11 and Theorem 14.21 together with (|1.4p we observe that 



(p e -p)\<u\ p = / ( Pe _^)(| ti |P_ c )+ / ( Pe -p)c 
n Jf2 Jo 

(5 ' 8) <C e ||V u ||^ (n) + Ce|| c || lv i. I , ( n) 

< Ce(||V U || p p(0) + |c|). 

Observe that u G W 1,p (O) ®I={h£ W 1,p (fi) : ti = c on dQ with c G R}. If 
u = v + c G WQ ,p (fi) © R, by the Sobolev trace inequality it follows that 

(5-9) \\<w^o) * CT\\u\\ p LP{gn} = c T \c\ p \dn\ p . 

Moreover, 

(5.io) h\\ p wl , m > \\Vu\\ p LP{ny 
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Then, by (fO]) and (|5TT0|) it follows that 

(5.11) HVu||£ p(n) + |c| < C\\u\\ 

From (JSU), (pUT]) and Lemma[0]it follows that if u S H^ 1,p (0) © E then 
(p e - p)M p < C||«||* (P1 „ (n) < C(G(«) + F(u, V)) 

where V is an arbitrary weight. 

Taking into account this remark, the proof in the non-flux case is analogous to 
the Neumann one. Observe that for this problem, we only need the boundary <9f2 
to be Lipschitz since we have used Theorem 14. II 

Robin: The extra term (3H(u) is irrelevant in the proof. This case is completely 
analogous to the Neumann and it follows by means of (|1.4p . (|2.ip . Theorem 14. 2 [ 
Lemma 14.31 and Lemma 14.41 



Eigenvalue depending boundary condition: From Theorem 14.21 Lemma 
and (|1.4I) we have that for every u £ W 1,P (Q) 

F{u, P) + H{u) < l |M| w -i. P(0) 



F(u, Pe ) + H{u) - F(u, p E ) + H{u) 

_G(u) + F{u,V) 



< 1 + Ce 



F(u,p) + H(u) 



and 



G(u) + F(u, V £ ) < G(u) + F(u, V) ^ Hk^n) 



H{u) + F{u,p) ~ H(u) + F{u,p) H(u) + F{u,p) 

G{u) + F{u,V) 



< (1 + Ce) 



#(u) + F(u,j5) 



Observe that (|1.4[) and (|2.ip allow us to bound in the opposite sense, that is, for 
each u e W^iQ), 

F(u, Pe ) + H(u) < G( M )+F( M ,T/ £ ) 
F(u,p)+H(u) ~ F(u,p e )+H{u)' 

G(u) + F(u,V) G(u)+F(u,V £ ) 
H(u) + F{u, p E ) ~ 1 + j ff(u) + p e ) ' 
Having changed this slight detail, the proof is analogous to the Neumann case. 

5.2. Proof of Theorem 11.21 The proof of the Stcklov case is simpler due to that 
in this case there is only an oscillating weight in the equation S e (d). 

Let 5 > and let G$ C W 1,p (£l) be a compact, symmetric set of genus k such 
that 

(5 . 12) Xt=sap m±f^n +0{S) . 

Being G\ admissible in the variational characterization of Af , we get 

(OS, K < S u P ofio+wa. 

ueQk H{u) 
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For every function u £ Gg C W 1,p (Sl) we can apply Theorem 14.21 obtaining 
G(u) + F(u,V s ) G(u) + F(u,V) j| 

(5 ' 14) #W #W ■ 

Now, for each u £ Gt we can apply Lemma l4~3l to bound 



\\\u\ p \\w^(n) G(u) + F(u,V) 
H(u) ~ C H(u) 
(5.15) < c sup G(,) + F(t,,V) 

= c(A fe + 0{S)). 

Then, combining (|5TT3|) . (|5~T4|) and (f57L5|) we find that 

A| < X k + 0(5) + Ce(X k + 0(6)). 

In a similar way, interchanging the roles of Xk and A|, we can obtain the opposite 
inequality Letting 5 — > we arrive at 

|A% — Afc| <Cemax{A fe ,Af}. 
By Lemma 13.21 we obtain that 

|A|-A*| <Cek^ 

and the proof is complete. 
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